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ABSTRACT 


This report deals with the development of a one dimensional steady state 
stratospheric aerosol model and the subsequent perturbations caused by 
Including the expected space shuttle particulate effluents In the model. Two 
approaches to the basic modeling effort have been made: In one, enough 

simplifying assumptions were Introduced so that a more or less exact solution 
to the descriptive equations could be obtained; In the other approach very 
few simplifications were made and a computer technique was used to solve the 
equations. The most complex form of the model contains the effects of 
sedimentation, diffusion, particle growth and coagulation. Results of the 
perturbation calculations show that there will probably be an immeasurably 
small Increase In the stratospheric aerosol concentration for particles 
larger than about 0.15 urn radius. The increase In very small particulates 
(greater than .01 urn diameter) is potentially large but cannot be adequately 
evaluated until the true natural background cf these small particles is 
determined. 


itratospheric Aerosol Model with Perturbations Induced by the Space Shuttle effluents 

J. M. Rosen, D. J. Hofmann, and S. P. Singh 

Department of Physics & Astronomy, University of Wyoming 82071 

I ntroduction 

It Is well known that a relative maximum exists on a global scale 
In the aerosol mixing ratio for particles having a diameter >0.3um approximately 
10 km above the tropopause. The particles are composed of H^SO^ (for the last 
decade at least) with many of thein containing smaller solid Inclusions. Following 
a large volcanic eruption the aerosol layer may experience a significant Increase 
In concentration; but during quiet periods of volcanic activity the layer In 
principle could approach a quasi steady state distribution. It Is the purpose 
of this paper to present a one dimensional model describing this quasi steady 
state situation and to Investigate perturbations of this state. 

The Model 

At the heart of the nx)del Is the assumption that a supersaturated layer of H^SO^ 
exists about 10 km above the tropopause. Although the model Itself Is not 
concerned with the chemistry of this layer's source. It could be formed from 

sulfur bearing gases such as SOo or CSO diffusing up through the tropopause and 

after a series of chemical reactions eventually forming H^SO^; or the source could 
be continual small volcanic eruptions with the required effective Injection 
altitude. Since H^SO^ has a very low vapor pressure at stratospheric 
temperatures and water vapor concentrations (Gmitro & Vermeulen, 1%3) even 
a very modest production rate of H 2 S 0 ^ could produce a large super saturat ion. 

The model further assun»es that this saturated vapor condenses on any 
particles that are present at a rate governed by the thernwl flux of H^SO^ 

molecules onto the particle.s' surface. Replenishmont of the particles comes 

from diffusion of tropospheric aerosol upward and as an option other 
sources can also be Included. The effect of coagulation Is also taken into 
account by the model . 


In this present work evaporation of the particles is neglected. This 
is more or less Justifiable because such a process would eventually saturate 
the entire stratosphere with H^SO^ vapor providing there were no sinks for 
the vapor. Under saturated conditions neither growtii nor evaporation would 
take place. However at high altitude there may be a sink for H^SO^ due to 
photodissociation by ultraviolet light. Under this circumstance the model 
presented here should not be considered entirely realistic above about 30 
to 35 km. On the other hand this is above the main region of interest and 
should have little affect on the major results presented. The inclusion of 
evaporation effects would increase to a considerable extent the complexity 
and uncertainties in the model. The basic reason for this is that the 
particles cannot evaporate to a snuller size than that of their original 
core. Since individual particle identities are lost in the diffusion 
process, the original core size is unknown. 

As noted above, the H^SO^ profile is treated as an adjustable parameter 
rather than deriving it from an appropriate chemical reaction model. 

Although this approach at first may seem unrealistic and questionable, the 
assumptions involved in present sulfur chemistry ntodels (including the 
nvjgnitude and type of sources) are simply too uncertain to produce a reliable 
HpSO^ vapor profile. Our interpretation of a reasonable profile is one which 
has a relative maximum in supersaturation near the observed aerosol nviximum, and 
a concentration consistant with the over all sulfur budget and transport propertie*; 
of the atmosphere. For this reason we have chosen to make the H^SO^ vapor 
profile a paranx?ter that can bo adjusted to obtain a good model fit to the 
observed aerosol profiles. It should be expected that similar results would 
be obtained from a more elaborate model containing the sulfur chemistry provided 
the HpSO^ vapor profiles in each case are the saim>. Such a comparison has boon 
made and will be discussed in a later section. Again it should be emphasized that 
in our opinion the present uncertainties involved in the sulfur chemistry make 


a nwre elaborate model (one containing sulfur chemistry) less fundamental than 
the approach presented here. 

Philoso phy o f one Dime nsional Models 

In general the atrwsphere should be treated as a 3 dimensional system 
and the applicability of one dimensional models Is open to considerable debate. 
Under present circumstances the limited amount of detail that can be Included 
In 3 dimensional models severely detracts from their credibility. Thus 
regardless of the number of dimensions used In a model, the applicability of 
the results will be open to a certain amount of justified criticism. In order 
to overcome this dilemma It seems reasonable to make a working assumption that 
one dimensional models can presently be used for a testing or a proving 
ground of new Ideas and concepts. If such Ideas show promise they will naturally 
(and eventually) evolve Into accurate multi dimensional models. Another 
useful working assumption Is that a good one dimensional model will capture 
some of the main general features of the constituents of interest. In the 
case of stratospheric aerosol this assumption is at least plausible because it 
Is known from direct measurement that the aerosol does in fact have recognizable 
features on a global scale. Obviously, in comparing one diniensional model 
predictions with actual field data, agreement between the oeneral character of 
the profiles is more important than an exact absolute quantitative agreemtot in 
a limited region. This view will be adopted here in comparing model predictions 
with typical field measurements. 

Basic Mode l Equ a U on 

The basic equation describing the time rate of change of the differential 
size distribution n(r,z) at altitude z is 

, IL -i Gnl+C 

3t 3z 3r 

where r is the radius, F is the particle flux (due to sedimentation and eddy 
diffusion). is the growth rate, C is the coagulation term and n is the 
particle number concentration. 


T h e SJ ny I i f 1 ed App roa c h 

Under certain conditions the above equation can be solved by s1m('1e 

analytical wthods. This Is done by usinq the equilibrium condition 

(.'>n/St*0), taking the eddy diffusion coefficient D to be constant 
2 

(bOOO cm /sec In most cases), letting C*0. using a A - 'unction for G, 
approximating the gravitational settling speed by a function that Is 
proportional to particle radius and Inversely proportional to air density and 
requiring an Isothennal stratosphere. Although some* of these simplifications 
may sewi quite crude It will later be seen, after comparison with solutions 
containing more realistic detail, that this simple approach captures the 
essential characteristics of the model both qualitatively and quantitively. 

The boundary conditions are specified by a source of single size "seed particles", 
corresponding to condensation nuclei (cn). at the lower bounda»y (tropopause) 
which diffuse »p to tlie 5- function q?*owth layer, a sink for these particles, 
for simplicity the specific gravity of the seed particles Is taken to be me 
same as that of H^SO^. Larger particles are formed at the growth layer and 
distribute themselves under the Influence of eddy diffusion, sed1i’x*ntat ion, and 
a sink at the tropopause. A schematic diagram sunxiurizlng these proc»*sses can 
b»* found In Appendix I. 

In this paper it Is assumt'd that every molecule of that strikes 

an aerosol particle will stick and liiriedlately 2 water vajKir molecules are 
taken on by the particle. This will result In about a 75* solution of 
for the stratospheric aerosol which Is In agreenx'iit with measurement (Rosen, 1R71). 
Growth due to co’llslons between M2S0^ nxilecules has been shown to be negligible 
in the stratosphere (Hamlll et al, 1977) and has been neglected. 

With tiv above simplifications the solution Is mathematically very 
similar to tliat obtained by simply requiring a iS- function source at the 
desired aerosol maximum, but with one Important addition: It provides tiie 

mechanism for generating a size distribution of stratospheric aerosol from 


the sincjle size seed particles diffusing up from the tropopause. 


The absolute concentration at the aerosol maximum is in part determined 
by the total number of H^SO^ molecules in the ^-function growth layer and 
the concentration of seed particles at the tropopause. In this model the 
flux of seed particles into the stratosphere is balanced by the flux of 
all the larger particles out of the st '-^tosphere. More detail concerning 
the derivation of the analytical form of the solution is presented in Appendix I. 

One minor and correctable shortcoming of this solution is tri? requirement 
of only one size particle at the tropopause when in reality an entire size 
distribution should exist. This fault can be alleviated by summing the 
solutions generated by a series of values for the seed particle sizes that 
approximate the desired size distribution at the tropopause. The result of such 
an approach will here be referred to as the quasi analytic solution and can be 
made as accurate as desired by increasing the number of points needed to 
approximate the size distribution at the tropopause. 

Resul ts 

A comparison of solutions using one size for the seed particles at the 
tropopause and using a complete size distribution at the tropopause is shown 
in figure 1. The insert shows the actual two integral size distributions used. 
Appropriate parameters have been chosen to reproduce the general character of 
the observed aerosol profiles (as shown in figure 4). The particles referred 
to as cn are actually the total number of particles present above a diameter 
of .01 nm and the profiles that refer to aerosol are the total number greater 
than 0.3 urn diameter or those generally associated with the stratospheric 
sulfate layer. The remaining profile is the ratio of the concentration of 
particles with diameter greater than 0.3 pm to the concentration of particles 
greater than 0.5 p diameter. As can be seen, there is barely a significant 
difference between the two examples and most of this can be attributed to the 
difference in average size of the two classes of seed particles. In the cn 
profile the particles associated with the smooth size distribution cannot 


diffuse to as high an attitude due to their larger average sl/e and the vonseguent 

greater Influence of sedlrientat Ion, In the aerosol profile, particles assoclateil 

□ 

with the sim^oth sire dlstrlhutlon are more numerous because they start out at a 
lartjer average sire, that Is, for a given anniunt of gri>wth nuMV of them will ivach 
0.3 kiin dlami'ter than the other class of particles. It can be generally concludeil 
that the results are not very dependent on the exact seed ivirticle sire distribution 
used at the t»*opopause (lowt*r boundary) and that using a sire distribution rather 
than a single sire Is In most cases not worth tlu' considerable extra effort and 
com|ilex1tv. 

The complete sire distribution at the aerosol nMxImum generated by using 
single sire .01 micron dianvter seed particles at the tropopause Is shown 
In figure ?. For ciwiparlson, results of experimental measurements are also 
shown In the figure. This exper1rx*ntal data has been critically dlscussini In 
somi‘ detail by Harris and Kosen (1?76). It Is apparent that the agreetm'nt between 
the sln^'le nx'drl aiKl the data Is surprisingly good and wi'uld tend to support the 
cn'dlblllty of the n»odel1ng approach. In addition, the mass flux of sulfur lunnied 
to sustain the lUSO, vapor laver is abiHit 3.3 x 10^ tons'vr. and i>, reasonably 

i. H 

consistent with otht*r estimates (4 to 17 x 10*^ tons/vrl as discussed by Trutren 

4 

(l‘T7b). A sulfur flux of 5 x 10 tons yr was used In Crutren's l*S0 nxHlel of 
St ra tospheric aerosol s . 
i'erturbat Ions 

Assuming that this simple approach to a solution captuivs tiv essence of 
tin' real stratospheric aerosol, it is of sow Interest to examini' the conseguence 
of several types of perturbations. Figure 3 shows the result of a perturbation 
causiHl by a nx'teoritic dust soutve and a space shuttle Al^O^ particle souixe. 

Only the cn profile Is slniwn because the change in the .1 » 0.3 urn aerosol and 
ratio profiles was found to be only of the order of lOt. These perturbations 
were calculated under the conditions of a constant number of I'^SO, molecules in 

4 , 4 


the growth layer. It would also be of interest to perform the calculations with 
a constant source strength of the H^SO^ vapor In the growth layer. Under this 
latter condition the resulting perturbations due to growth would not be as large 
because the number of HpSO^ molecules the growth layer would not be as great. 
Thus the results shown In figure 3 were calculated under assumptions which would 
produce the largest effect. 

Mathematically the meteorltic dust source Is treated similarly as the 
tropospheric source of seed particles. The upper boundary Is taken high 
enough so that only sedimentation and diffusion are Important processes. 

The meteorltic particle concentration Is chosen to be consiste't with a 
conservative meteorltic flux (In this case, 10^ tons/yr with an average 
radius of .04 microns and an average density of 2 gm/cm^). Thus the growth 
layer acts as a sink for both the tropospheric seed particles as well as 
meteorltic particles and the net flux from both of these sources Is the 
relevant quantity to be used In the balance equations discussed In Appendix I. 

The space shuttle oerturbatlon Is dealt with as a superposition of 
solutions. Each Individual solution Is that for a point source at the 
desired altitude Increment. Thus the net profile Is a sum of a large 
number of exact solutions. The absolute values of the Injection rates 
are the same as those used by Hofmann et al (1975). 

For reference, the expected cn profile for the space shuttle only 
(no growth layer, tropospheric source or meteorltic source) Is also shown 
In figure 3. It Is Interesting to note that the space shuttle perturbed 
profile Is practically the sum of the profile obtained for the spate 
shuttle alone and the natural cn profile predicted by the model. TMs 
result along with the previously mentioned fact that the d > 0.3 um aerosol 
profile is practically unchanged under the Illustrated perturbations suggests 
that the equilibrium distribution of the space shuttle particulate effluents Is 
practically Independent of the aerosol growth modi ’ here. Thus It would 

appear that simple calculations neglecting th® g mlcs of the natural 

aerosol may be reasonably valid. 


It should also he noted that the basic parameters used In the model 
to generate f Inure .1 are not exactly the same as those used to qcnerate flqurt* 1. 
The principle difference is a reduction in the t .iicentration of particles at 
the tropopause by a factor of 0.4 in the latter figure. This change tends to 
bring the absolute aerosol concentration at the maximum into better an»*eenx»nt 
with the measurements (shown in figure 4) but the general character of the 
profiles remains unchanged. 

The result of this simple perturbation study indicates that the space 

shuttle could change the high altitude cn concentration by a large factor if 

meteoritic sources are not important. However even a small metec itic source 

could change this conclusion. In addition, higher values of the coefficient of 
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eddy uiffusivity than those used here ^5000 cm /sec) would also reduce the effect 
of the space shuttle at high altitude. At present, good experimental measureixMits 
are needed above km to determine typical cn profiles, and until these 
measurements are made it will not be possible to determine the extent to which 
the space shuttle activities will increase the cn concentration in the upper 
stratosphere. According to the prediction of this simple model the effect of Uh’ 
space shuttle on the d ^ 0.3 pm aerosol profile appears to be quite negligible. 
Computer Solution Approach 

Of all the approximations used to obtain a quasi analytic solution, only 
two may be questioned as not being realistic: the assumption of a 5-function 

growth layer and the neglect of coagulation. Fror: the work of Hofmann et al (1976) 
it is deal that a fairly narrow source region would be required to successfullv 
explain the observed stratospheric aerosol layer. The use of a 5-function to 
describe the source region is therefore not entirely unrealistic. Thus the 
principle short coming of the quasi analytic solution is the complete neglect of 
coagulation. In what follows, a more general and more complete solution to the 
basic equation will be obtained by employing computer methods. 
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A finite difference *4thod has been used to obtain a computer solution to 
the basic equation. The altitude grid size for most cases was taken to be 1 km 
and the radius grid points differed by factors of essentially 2^^^. The method 
of solution Involved starting with Initial conditions and letting the 
differential equation develop In time until a steady state solution was reached. 

In most cases, time Intervals of one day were used. Errors that develop In this 
method due to finite grid size were Investiqated by decreasing the grid sizes and 
comparing the resulting successive steady state solutions. 

The eddy diffusion profile used Is shown In figure 4. The upper portion 

Is similar to that suggested by Hunten (Johnston et al , 1976) but the lower 

portion has been modified to that of Chang (Johnston et al, 1976) which we feel 

Is consistent with a tropopause at 10 km. Recently other modelers have adopted 

a profile similar to the one used here (Luther 1977). We have observed that 

the resulting particle profiles are not very sensitive to the exact nature of 
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the eddy diffusion profile and the use of a constant value of 5000 cm /sec would 
not change the character of the solutions which will be Illustrated. 

The gravitational settling velocities have been adopted from Kasten (1968) 
by fitting his tabulated data to a curve that Is proportional to the particle 
radius and specific gravity and Inversely proportional to the ambient air density. 
A particle specific gravity has been used that Is consistent with a composition 
of about 75‘ H^SO^ and 25* H^O. The vapor pressure profile of H2S0^ was taken 
to be a gausslan function centered at 20 km and several kilometers wide. A 
peak H2S0^ concentration of 2.75 x lO’^^ mm Hg was chosen, a value In considerable 
excess of 100‘. saturation. 

The cumulative size distribution at the tropopause was taken to be 
Inversely proportional to the 3.5 power of the radius and specified by an 
absolute value that was consistent with field measurements of the aerosol. 

A smooth lower cut off In the size distribution near .05 lim radius was 



found necessary In order to obtain aqreoment In the absolute values of tin* 
calculated and measured aerosol profiles. A$ will be seen this cut off leads 
to using smaller values of the cn concentration at the tropopausc than are 
actually observed. The significance of this problem will be discussed In a 
later section. 

The upper boundary was high enough so that the flux of particles across It 
could be taken as zero (If a meteorltic source was Included then the flux was 
derived from the corresponding concentration at the level of the upper boundary). 

Due to limited computer facilities an approximation to the treatment 
of coagulation was developed. A constant value of the coagulation coefficient 
K was used (36 x 10*'® cm^/sec) that was obtained by averaging over typical 
size dfstrlbutlons encountered. It should be pointed out that an Inconsistancy 
of a factor of two <n the definition of K has developed In the literature. We 
use the notatio' ..nu formulation of Walter (1973). The values reported b\ Fuchs 
(196‘*) for Instance are about a factor of two lower which can be attributed 
to a different definition of K. 

An expression for the time rate of change of the size distribution due 
to coagulation (which Is required In the basic model equation) has been given 
bv U. Iter (1973). It was necessary to develop an approximation to this 
expression consistent with the radius grid size. Since the details of tlie 
required approximat ion are cumbersome and tedious, they will not be presented 
here. In its place, a discussion of the overall model consistency checks will 
be nude. 

S olution Checks 

Since there are many opportunities for serious errors to develop in the 
computer solution (i.e. those due to approximat Ions , cumulative errors and 
programming mistakes) it is essential to have some independent means of 
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verifying the over-all results. Conservation of mass was checked by conniarlnq 
the mass loss of vapor to the net mass loss of the particles diffusing 

out of the stratosphere (It Is necessary to make a correction to the latter 
mass flux due to Its partial water content as previously described). This 
test Is very sensitive to the accuracy of the treatment of coagulation. In 
the limit of no coagulation the solution can be checked for particle conservation 
because the flux of seed particles Into the stratosphere must be balanced by 
the flux of larger particles out of the stratosphere when equilibrium Is established. 
Also In this limit It Is possible to judiciously choose the parameters so that a 
direct comparison between the computer solution and the accurate quasi analytic 
solution can be made. However, In this case good altitude resolution Is needed to 
reasonably approximate the 6-functlon H^SO^ vapor profile with a narrow qausslan 
curve required by the computer method. Results of these tests Indicate that thc» 
overall accuracy of our computer solution Is about lOT.. 

Results 

Typical computer generated profiles are shown In figure 4 along with the 
range of actual measutro values for a tropopaure near 10 km. The experimental 
data was taken from Hofmann et al (1975), Rosen et al (1975), Pinnick et al 
(1976), and Rosen and Hofmann (1977). 

The size distribution at the aerosol maximum associated with the profiles 
shevn In figure 5 is presented In figure 2 and Is very similar to that of the 
quasi analytic solution. The mass flux of sulfur required to sustain the 
vapor layer, used In calculating the profiles of figure 5, is about 8 x 10 ton*, 
per year and Is In reasonable agreement with that obtained from the simple iix)del 
1i lustrated In figure 1 . 

Even though the computer model contains the effects of coagulation and 
the quasi analytic model does not, the two are In reasonable agreement. This 
result can be attributed to the low particle concentration used at the tropopause 
which subsequently yields only slight coagulation affects. The Influence of 
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at high aUltuite. On thn other hand evaporation would art as a s^.nik for largo 
particles which explains the lowi'r values of aerosol concentration In the Ainos 
model at high altitude. 

It should h«* |»o1nted out that a rigorous comparison between thi‘ twti 
nwdols Is not strictly possible binause they do not require Identical 
Input paraitH'ters. Son»» Judgemt^nt Is required as to what exactly would 
constitute a mutually consistent set of paraim*ters. Thus the results 
shown In figure 6 are not entirely objective. Howi'ver, within the range of 
jK)ss1ble mutually consistent paranx’ters. It Is our experience that the two 
nwdels prtxluce pixiflles that are very similar In character. 

C onclus ion 

Although many asfH*cts of tht> stratospheric aerosol are consistent with 
this nxulel.the predicted cn profile Is not In good agreiMnent with observation. 

It Is possible to change the Input parameters In such a way so as to bring the 
cn profile into better agreement with rxMSureim'nt s (see figute b for example) 
but this would »'«*sult In other iinacceptable profiles (I.e. the d 0.3 pm aerosol 
profile). We did not find a reasonable set of paraim'ters that would bring the 
itxulel intc essentially iwnplefe agreement with the measurements. 

The disagreeim*nl between the predicted and im'asured cn profiles at high 
altitude (figure 1 and ft) Is probably not serious. This discrefxincy could be 
corrected by using large?* values of eddy diffusion at high altitude (as was 
done ii? figure 6) or adding a iix*teoritic source, the effects of wl\ich a?v 
illustrated in figute 3. Itv' problem with the cn profile near the 
tropopause si'ems to be im>re basic: the nKHlel simply does not recjuire as nuny 

cn as are measu.vd. 






In conclusion wt> feci that the model presented In this paper Is capable 
of describing most of the features of the stratospht»r1c aerosol but certain 
shortcomings Indicate that the model Is not entirely correct or complete. 

If these shortcomings are overlooked aitd the model Is used to predict the 
change In particulate content of the stratosphere due to the space shuttle 
Al^Oj effluents, then only a small and probably Immeasurable change In the 
concentration of particles larger than 0.3 urn diameter would be expected to 
occur. Although a potentially significant amojnt of smaller particles 
(d ^ .01 urn) will be added to the stratosphere above 30 km It Is not 
possible to determine If a measurable Increase above natural levels will 
occur simply beacuse tin? natural background In this size range and altitude 
Is unknown. 
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FIGURE CAPTIONS 


Figure 1. A comparison of profiles using two types of size distributions. 

The insert shows the actual cumulative size distributions employed. 

The cn profile refers to particles larger than 0.01 urn diameter and 
the aerosol profile refers to particles larger than 0.3 nm diameter. 

The ratio profile refers to the ratio of the concentration of particles 
with diameter greater than 0.3 urn to the concentration of particles 
with diameter greater than 0.5 urn. 

Figure 2. A comparison of tne size distributions at the aerosol maximum 
derived from the quasi-analytic solution (thick solid line) and the 
computer solution (heavy dotted line) with the actual measurements. 

The difference in the two types of solutions at small radii is due to 
the different choice of seed particle concentration assumed at the 
tropopause. See table for key to references. 

Figure 3. The cn profile obtained from tne quasi analytic solution using 
various assumptions: TS ■ tropospheric source (of seed particles), 

G * (S- function growth layer at 20 km present; MS ■ mcteoritic source 
(of seed particles) and SS » space shuttle source (of seed particles). 

Figure 4. The eddy diffusivity profile used in the computer solution compared 
with a profile suggested by Hunten. 

Figure 5. A typical computer solution (solid lines) compared to the range 
of actual measurements (dashed lines). The cn profile refers to 
particles larger than .01 um diameter and the aerosol profile refers 
to particles larger than 0.3 pm diameter. The aerosol ratio refers 
to the ratio of the concentration of particles with diameter greater 
than 0.3 nm to the concentration of particles with diameter greater 
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than 0.5 ;mi. 


Figure 6. A compArlson of the Ames Model (dashed lines) and the Wyoming 
Model (solid lines) using similar Input paraiiieters. See text for 
explanations of differences. The cn profile refers to particles 
larger than .01 urn diameter and the aerosol profile refers to particles 
larger than 0.3 urn diameter. The ratio profile refers to ratio of the 
concentration of particles with diameter greater than 0.3 urn to the 
concentration of particles with diameter greater than 0.5 urn. 




TABLE 

Figure 2 Reference Code 


Code 

Reference 

Type of Measurement /Comments 

B 

Bigg (1975) 

Impactor 

BFT 

Brownlee, Ferry, & Tomendl (1976) 

j Impactor for large particles 

1 

F 1 

Friend (1966) 

j Impactor 

FL 

Ferry & Lem (1974) 

Impactor 

I 

Ivlev (1976) 

Impactor/Data from Aug. 1975 

JCM 

Junge, Chagnon & Manson (1961) 

Impactor/Data from 26 Aug. 1958 
(318.4 km 

M 

Mossop (1964) 

Impactor/Data from 7 April 1964 

MD 

Miranda & Dulchinos (1975) 

Photoelectric Particle Counter 

MOM 

Miranda, Dulchinos A Miranda (1973) 

Photoelectric Particle Counter 

W 

Wyoming Results (Pinnick et al 1976) 

Photoelectric Particle Counter 
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Appendix I 


The niethod of obtaining a solution to the simplified sedimentation • 
diffusion • growth equation, as discussed In the main text, can be easily 
understood by referring to figure A-1. The 6-funct1on growth layer acts 
as a sink for the seed particles that originate at the tropopause. The 
flux of these particles to the growth layer Is therefore easily found to 
be (see example 1, Appendix II): 

♦in ■ P(Z(,) ur e*'^s/^m (l-e*'"s^^m)‘^ 

where 

p (Zp) ■ density of air at the tropopause. 

u » constant In the sedimentation velocity equation: v ■ ure^^*^ 
r » radius of particles 

Tj “ mixing ratio at tropopause of beed particles. 
z_ 

1/q « (e -1) UH/D 

m 

H ■ scale height of atmosphere. 

D ■ eddy diffusivlty (constant) 

■ height of 6-funct1on growth layer above tropopause. 
r^ ■ radius of seed particles. 

Consider the particles in the growth layer itself. In each size interval 
there will be a gain In concentration due tt smaller particles growing larger, 
j loss due to particles growing out of the size range and a loss due to sedi- 
mentation and diffusion to the tropopause. The growth can be treated as a 
"flux" of particles through the size distribution given by Gn where G Is the 

growth rate (In radius units/time) and n is the differential number concentration 
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Figure A-1. A schematic diagram of the stratospheric aerosol model. 
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(cm'’’ liin" ). The "divergence" of this "flux" must therefore be the loss by | 

sedimentation and diffusion to the troposphere. From the arguments given ' 

In Appendix II (see example 4 ) this loss term Is found to be: 

‘^'^out * ^0^"" 

where Is the equilibrium differential mixing ratio In the growth layer. 

On the other hand the loss term from the "divergence" of the "flux" of particles 
through the size distribution Is 

I 

i 



Under present circumstances G Is due essentially to the thermal flux of H^SO^ : 

onto the particles. This is given by G ■ '4 N v^^^ where N ■ the concentra- f 

tion of H2S0^ molecules, v^^ ■ the effective velocity of H^SO^ molecules and | 

Vg Is the effective volume of 1 H2S0^ and 2 H2O molecules. The units of | 

are number of particles per unit volume In size interval dr being lost per | 

second. Thus the losses of a given size particle can be found by Integrating 
over the altitude range of the layer Itself. Since the layer here is a 
6-function the result Is 

*g ■ ''th "e fr 

where A = total number of H2S0^ molecules in the layer. 

The total flux of particles out of the stratosphere can be found by Integrating 

; 

the expression for d 4 >Qy^ over all size ranges; this quantity will be referred 
“ “ ♦out- 

d 

By equating 4 *^^^ to 4 >g a simple differential equation Is obtained which 
yields a relative size distribution in the growth layer (i.e., at z^^). 


By equating to the absolute size distribution Is obtained (that Is, 
Y(z^,r)). The value of y at other altitudes Is obtained by simply treating 
the growth layer as a source (with a fixed mixing ratio y(z„,*>^)) and using 
sedimentation-diffusion equations to determine the spreading (see Appendix II). 
The results of this calculation are as follows: 

At the aerosol maximum the mixing ratio Is given by 


Hz„.r) 


fr. 


[ar ♦ 


2 3 

br^ ♦ cr”* + 


where 


-z /H -1 

a - 4D/(v^ Vg AH) [1-e ] 


u _ a ^ H r 
b • 4 - -y [e 


Zm/H 

- 1 ] 


c » 


0.4b' 


-4r.b/a -4r b/a -1 , lo* ^ , 

f«4r5be ^ (1-e ^ V 


r^“ mixing ratio of seed particles at tropopause. 


r^» radius of seed particles. 


*Approximate value for range of values of r^, a and b expected in the 
stratosphere. 

Only a few terms are needed in the exponent describing y(z„,*'^) and in 
fact Including only the first term is generally adequate for the size ranges 
in which the concentrations are high enough to easily measure. This result 
shows that a stratospheric aerosol generated by a growth process is probably 
better described by an exponential size distribution than by a power law size 
distribution. 


A sample calculation using parameters Involved In generating the profiles 
shown In figure 1 Is given below. 

v^. ■ i/SKT/mn ■ 2.15 X 10^ cm/sec (T ■ 216° K and m ■ mass of 

H^SO^ molecule) 

Vg ■ 1.3 X lO’^^ cm^ (approx. 75* H2S0^ 25* H 2 O) 

A ■ 1.47 X 10^^ H-SO* moluecules/cm^ (Compare this value to 3 x 10^^ 

^ ^ used In the Computer Model 

results In figure 5) 

Vj ■ 8.9 X 10 “^ r e^^^» (r In urn) 

H ■ 6.3 km 

D » 5000 cm^/sec 

■ 10 km 
m 

r^ ■ .01 urn 

3 

p (z) ■ 4.127 X 10' e (z ■ distance above tropopause) 

■ 1.21 X 10^/gm (or 50/cm^) 

With these parameters the following results were obtained: 


1 

function 

value 

a 

9.71 

um ^ 

1 

b 

10.59 

um'^ 


4.62 

um'^ 

f 

10.63 


N(.15) at z^ 

1.32 

cm'^ 

N(.15)/N(.25) 

5.0 



The complete profiles are shown In figure 1. N(.15) and N(.25) refer to 


the concentration of particles larger than .15 um radius and .25 jim radius 
respectively. 



Appendix II 


Library of Profiles for Simple Diffusion- 
Sedimentation Equilibrium. 

It Is often quite Informative to refer to simple one-dimensional eddy 
diffusion-sedimentation profiles of aerosol In order to gain an approximate 
understanding of the spatial distribution of this contituent. In addition, 
complex profiles can often be viewed as a combination of several simple 
profiles. 

Under equilibrium conditions the flux of particles of radii between r 
and r-t^dr and at altitude z Is given by 

♦(z.r) • -p(z) v(z.r) i(z,r) - p(z) D(z) 

S ■ ■ Iz ^ ® 

where 

p(z) ■ air density 
v(z,r) ■ fall velocity 
y{z*r) ■ mixing ratio 
D(z) » eddy diffusivlty. 

These equations can be solved exactly for the following conditions which 
are fairly realistic. 

o(z) = Pq (H ■ atmospheric scale height) 

v(z,r) ■ ure^^^ (u * constant) 


D(z) * constant 


Under these conditions a general solution Is: 


Y(z.r) ■ f(z.r) ♦ Yj, 

where and y^^ are constants determined by the boundary conditions and 

f(z,r) - exp[-(e^/“-l) Hur/D] s 

when the lower boundary Is at Z"0* 
and 

f(z.r) . exp[-(e'/"-e'"/") Hur/D 5 e"-/9’(z) 

when the lower boundary Is at > 0. 

By definition 

l/g(z) = uH/D 

l/g'(z) = uH/D 

Some useful ' nmples of solutions follow. 

1. Source at z*0 and sink at z*z„ > 0. 

m 

y(z.r) ■ y(o.r) 


* ■ constant ■ o(ZQ)ur y (o,r) 



y ( 0, r ) 

MIXING RATIO (y ) 


2. Source at z*0 and sink at z*«. 


Y(i.r) ■ 



MIXING RATIO (/) 


3. Source at z»z„ and sink at z«®. 
m 

Y(z,r) • Y(z^.r) ; y(z^*r) given 


♦ ■ 0 



MIXING RATIO (y) 



i— i-J>— 










Source at z, and sink at z»z„. (z, > z„ >0) 

c n ctn 

y( 2 .r) ■ y(z 2 ,r) (1-e"'"'^^ 

♦ ■ - y(z2»>^) ^^2^)“^ 

Zg L — s 


SOURCE AT Z*Z2 


SINK ATZ'Zm 


MIXING RATIO ( /) 


Source at z ■ +«> and sink at z„. 

m 

Y(Ztf') ■ y(®.<^) (1 - yJ^.r) given 


♦ ■ constant ■ -y(®,r) PjjUr 



SOURCE AT Z = 00 


SINK ATZ = Z, 


r ( ®,r ) 

MIXING RATIO (r) 





